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Abstract 

Critical remarks in respect to the Generalized Relativistic Effective Core 
Potential in the recent article of S. A. Wildman at al. are discussed and 
shown to be incorrect. 



S. A. Wildman, G. A. DiLabio, and P. A. Christiansen's article |l] deals with a problem 
of the Relativistic Effective Core Potential (RECP) generation for precise calculations. In 
this article, new RECPs for the Tl, Pb, Bi, At, and Rn atoms which are more accurate 
than the previous RECPs |2| are constructed and tested in molecular calculations. However, 
some statements made in the article are incorrect. 

As it follows from article JTJ, the problems with the old RECPs |2| "... are due to the 
inappropriate partitioning of the spinors used to generate the f$/2 and f 7 / 2 potentials" that, 
in its turn, is connected with the restriction on the number (n < 2) of inflections for pseu- 
dospinors in the RECP generation procedure developed in |J. However, this does not mean 
that "Given the comparable approaches in terms of orbital/spinor partitioning, one would 
expect to see very similar behavior and problems with the potentials developed by other 
groups 16 . . ." (i.e., Ref. |§] here). Generalized RECP (GRECP) generation procedure f§] is 
based on procedure || but has essential distinctions. Some of them are requirement of a 
smooth shape for the potential and absence of limitations on the number of inflections for 
pseudospinors ||. We did not concentrate on such technical details in our papers, because 
they are not so important as the principal questions concerning the optimal form of the 
RECP operator. One can see from comparison of Fig. |I| in this comment with Fig. 1(B) in 
Ref. however, that the new 5/5/2 pseudospinor of S. A. Wildman et al., in contrast to the 
old one, is rather close to the 5/ 5 / 2 pseudospinor constructed within the GRECP generation 
procedure [||]. 

Some comments should be made on the authors' statement that in paper [Q "... it has 
occasionally been suggested that the small components are neglected in this procedure . . ." 
(i.e. in the GRECP generation procedure). More correctly (that is evident from our papers, 
see e.g. 0), we neglect contribution from the small components in the valence region to the 
two-electron integrals and the electronic density, and only the large components of Dirac 
spinors are used in order to construct the pseudospinors. In paper M, it was also specified 
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that we applied the two- component technique of Lee et al. || where the authors, in partic- 
ular, wrote at the end of page 5862: "Since small components become even less important 
in the outer region, it may be reasonable to assume that small components can be ne- 
glected in calculations that emphasize the description of valence electrons". Nevertheless, 
it is clear from this paper that the authors understand that the small components can not be 
neglected in Dirac-like equations. A few lines above in the commented article, S. A. Wild- 
man et al. write: "The two-component pseudospinor, is then effectively reinserted into 
the radial DHF equation and the equation is inverted to obtain the localized relativistic 
effective potential, Uj^ EF '(r)" . In reality, one should use "nonrelativistic-type" HF equations 
(i.e. containing the nonrelativistic kinetic energy operator) in the jj-coupling scheme which 
are augmented with the relativistic j-dependent potentials as it is made in ||. 
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FIGURES 




FIG. 1. The radial parts of the 5/ 5 / 2 pseudospinor (dashed line) and the large component of 
the 5/5/2 spinor (solid line) for the Pb atom from work [Q. 
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